Genetic networks underlying salinity tolerance in wheat uncovered with genome-wide analyses and selective sweeps.

KEY MESSAGE: A genetic framework underpinning salinity tolerance at reproductive stage was revealed by genome-wide SNP markers and major adaptability genes in synthetic-derived wheats, and trait-associated loci were used to predict phenotypes. Using wild relatives of crops to identify genes related to improved productivity and resilience to climate extremes is a prioritized area of crop genetic improvement. High salinity is a widespread crop production constraint, and development of salt-tolerant cultivars is a sustainable solution. We evaluated a panel of 294 wheat accessions comprising synthetic-derived wheat lines (SYN-DERs) and modern bread wheat advanced lines under control and high salinity conditions at two locations. The GWAS analysis revealed a quantitative genetic framework of more than 200 loci with minor effect underlying salinity tolerance at reproductive stage. The significant trait-associated SNPs were used to predict phenotypes using a GBLUP model, and the prediction accuracy (r2) ranged between 0.57 and 0.74. The r2 values for flag leaf weight, days to flowering, biomass, and number of spikes per plant were all above 0.70, validating the phenotypic effects of the loci discovered in this study. Furthermore, the germplasm sets were compared to identify selection sweeps associated with salt tolerance loci in SYN-DERs. Six loci associated with salinity tolerance were found to be differentially selected in the SYN-DERs (12.4 Mb on chromosome (chr)1B, 7.1 Mb on chr2A, 11.2 Mb on chr2D, 200 Mb on chr3D, 600 Mb on chr6B, and 700.9 Mb on chr7B). A total of 228 reported markers and genes, including 17 well-characterized genes, were uncovered using GWAS and EigenGWAS. A linkage disequilibrium (LD) block on chr5A, including the Vrn-A1 gene at 575 Mb and its homeologs on chr5D, were strongly associated with multiple yield-related traits and flowering time under salinity stress conditions. The diversity panel was screened with more than 68 kompetitive allele-specific PCR (KASP) markers of functional genes in wheat, and the pleiotropic effects of superior alleles of Rht-1, TaGASR-A1, and TaCwi-A1 were revealed under salinity stress. To effectively utilize the extensive genetic information obtained from the GWAS analysis, a genetic interaction network was constructed to reveal correlations among the investigated traits. The genetic network data combined with GWAS, selective sweeps, and the functional gene survey provided a quantitative genetic framework for identifying differentially retained loci associated with salinity tolerance in wheat.

Population-dependent reproducible deviation from natural bread wheat genome in synthetic hexaploid wheat.

Sequence elimination is one of the main mechanisms that increases the divergence among homoeologous chromosomes after allopolyploidization to enhance the stability of recently established lineages, but it can cause a loss of some economically important genes. Synthetic hexaploid wheat (SHW) is an important source of genetic variation to the natural hexaploid wheat (NHW) genepool that has low genetic diversity. Here, we investigated the change between SHW and NHW genomes by utilizing a large germplasm set of primary synthetics and synthetic derivatives. Reproducible segment elimination (RSE) was declared if a large chromosomal chunk (>5 cM) produced no aligned reads in more than five SHWs. RSE in five genomic regions was the major source of variation between SHW and NHW. One RSE eliminated almost the complete short arm of chromosome 1B, which contains major genes for flour quality, disease resistance and different enzymes. The occurrence of RSE was highly dependent on the choice of diploid and tetraploid parental lines, their ancestral subpopulation and admixture, e.g. SHWs derived from Triticum dicoccon or from one of two Aegilops tauschii subpopulations were almost free of RSE, while highly admixed parents had higher RSE rates. The rate of RSE in synthetic derivatives was almost double that in primary synthetics. Genome-wide association analysis detected four loci with minor effects on the occurrence of RSE, indicating that both parental lines and genetic factors were affecting the occurrence of RSE. Therefore, pre-pre-breeding strategies should be applied before introducing SHW into pre-breeding programs to ensure genomic stability and avoid undesirable gene loss.

Genetic and transcriptional variations in NRAMP-2 and OPAQUE1 genes are associated with salt stress response in wheat.

KEY MESSAGE: SNP alleles on chromosomes 4BL and 6AL are associated with sensitivity to salt tolerance in wheat and upon validation can be exploited in the development of salt-tolerant wheat varieties. The dissection of the genetic and molecular components of salt stress response offers strong opportunities toward understanding and improving salt tolerance in crops. In this study, GWAS was employed to identify a total of 106 SNP loci (R2 = 0.12-63.44%) linked to salt stress response in wheat using leaf chlorophyll fluorescence, grain quality and shoot ionic (Na+ and K+ ions) attributes. Among them, 14 SNP loci individually conferred pleiotropic effects on multiple independent salinity tolerance traits including loci at 99.04 cM (R2 ≥ 14.7%) and 68.45 cM (R2 ≥ 4.10%) on chromosomes 6AL and 4BL, respectively, that influenced shoot Na+-uptake, shoot K+/Na+ ratio, and specific energy fluxes for absorption (ABS/RC) and dissipation (DIo/RC). Analysis of the open reading frame (ORF) containing the SNP markers revealed that they are orthologous to genes involved in photosynthesis and plant stress (salt) response. Further transcript abundance and qRT-PCR analyses indicated that the genes are mostly up-regulated in salt-tolerant and down-regulated in salt-sensitive wheat genotypes including NRAMP-2 and OPAQUE1 genes on 4BL and 6AL, respectively. Both genes showed highest differential expression between contrasting genotypes when expressions of all the genes within their genetic intervals were analyzed. Possible cis-acting regulatory elements and coding sequence variation that may be involved in salt stress response were also identified in both genes. This study identified genetic and molecular components of salt stress response that are associated with Na+-uptake, shoot Na+/K+ ratio, ABS/RC, DIo/RC, and grain quality traits and upon functional validation would facilitate the development of gene-specific markers that could be deployed to improve salinity tolerance in wheat.

Decomposing Additive Genetic Variance Revealed Novel Insights into Trait Evolution in Synthetic Hexaploid Wheat.

Whole genome duplication (WGD) is an evolutionary phenomenon, which causes significant changes to genomic structure and trait architecture. In recent years, a number of studies decomposed the additive genetic variance explained by different sets of variants. However, they investigated diploid populations only and none of the studies examined any polyploid organism. In this research, we extended the application of this approach to polyploids, to differentiate the additive variance explained by the three subgenomes and seven sets of homoeologous chromosomes in synthetic allohexaploid wheat (SHW) to gain a better understanding of trait evolution after WGD. Our SHW population was generated by crossing improved durum parents (Triticum turgidum; 2n = 4x = 28, AABB subgenomes) with the progenitor species Aegilops tauschii (syn Ae. squarrosa, T. tauschii; 2n = 2x = 14, DD subgenome). The population was phenotyped for 10 fungal/nematode resistance traits as well as two abiotic stresses. We showed that the wild D subgenome dominated the additive effect and this dominance affected the A more than the B subgenome. We provide evidence that this dominance was not inflated by population structure, relatedness among individuals or by longer linkage disequilibrium blocks observed in the D subgenome within the population used for this study. The cumulative size of the three homoeologs of the seven chromosomal groups showed a weak but significant positive correlation with their cumulative explained additive variance. Furthermore, an average of 69% for each chromosomal group's cumulative additive variance came from one homoeolog that had the highest explained variance within the group across all 12 traits. We hypothesize that structural and functional changes during diploidization may explain chromosomal group relations as allopolyploids keep balanced dosage for many genes. Our results contribute to a better understanding of trait evolution mechanisms in polyploidy, which will facilitate the effective utilization of wheat wild relatives in breeding.

Allelic variations and differential expressions detected at quantitative trait loci for salt stress tolerance in wheat.

The increasing salinization of agricultural lands is a threat to global wheat production. Understanding of the mechanistic basis of salt tolerance (ST) is essential for developing breeding and selection strategies that would allow for increased wheat production under saline conditions to meet the increasing global demand. We used a set that consists of 150 internationally derived winter and facultative wheat cultivars genotyped with a 90K SNP chip and phenotyped for ST across three growth stages and for ionic (leaf K+ and Na+ contents) traits to dissect the genetic architecture regulating ST in wheat. Genome-wide association mapping revealed 187 Single Nucleotide Polymorphism (SNPs) (R2 = 3.00-30.67%), representing 37 quantitative trait loci (QTL), significantly associated with the ST traits. Of these, four QTL on 1BS, 2AL, 2BS and 3AL were associated with ST across the three growth stages and with the ionic traits. Novel QTL were also detected on 1BS and 1DL. Candidate genes linked to these polymorphisms were uncovered, and expression analyses were performed and validated on them under saline and non-saline conditions using transcriptomics and qRT-PCR data. Expressed sequence comparisons in contrasting ST wheat genotypes identified several non-synonymous/missense mutation sites that are contributory to the ST trait variations, indicating the biological relevance of these polymorphisms that can be exploited in breeding for ST in wheat.

Genome-wide association study for agronomic and physiological traits in spring wheat evaluated in a range of heat prone environments.

KEY MESSAGE: We identified 27 stable loci associated with agronomic traits in spring wheat using genome-wide association analysis, some of which confirmed previously reported studies. GWAS peaks identified in regions where no QTL for grain yield per se has been mapped to date, provide new opportunities for gene discovery and creation of new cultivars with desirable alleles for improving yield and yield stability in wheat. We undertook large-scale genetic analysis to determine marker-trait associations (MTAs) underlying agronomic and physiological performance in spring wheat using genome-wide association studies (GWAS). Field trials were conducted at seven sites in three countries (Sudan, Egypt, and Syria) over 2-3 years in each country. Twenty-five agronomic and physiological traits were measured on 188 wheat genotypes. After correcting for population structure and relatedness, a total of 245 MTAs distributed over 66 loci were associated with agronomic traits in individual and mean performance across environments respectively; some of which confirmed previously reported loci. Of these, 27 loci were significantly associated with days to heading, thousand kernel weight, grain yield, spike length, and leaf rolling for mean performance across environments. Despite strong QTL by environment interactions, eight of the loci on chromosomes 1A, 1D, 5A, 5D, 6B, 7A, and 7B had pleiotropic effects on days to heading and yield components (TKW, SM-2, and SNS). The winter-type alleles at the homoeologous VRN1 loci significantly increased days to heading and grain yield in optimal environments, but decreased grain yield in heat prone environments. Top 20 high-yielding genotypes, ranked by additive main effects and multiplicative interaction (AMMI), had low kinship relationship and possessed 4-5 favorable alleles for GY MTAs except two genotypes, Shadi-4 and Qafzah-11/Bashiq-1-2. This indicated different yield stability mechanisms due to potentially favorable rare alleles that are uncharacterized. Our results will enable wheat breeders to effectively introgress several desirable alleles into locally adapted germplasm in developing wheat varieties with high yield stability and enhanced heat tolerance.

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.).

Faba bean (Vicia faba L.) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely related crop species, lentil and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2-cM interval of Vf chromosome 2 which was highly colinear with Mt3. The obvious candidate gene from 78 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene controlling anthocyanin production in Medicago and resequencing of the Vf orthologue showed a putative causative deletion of the entire 5' end of the gene.

Global agricultural intensification during climate change: a role for genomics.

Agriculture is now facing the 'perfect storm' of climate change, increasing costs of fertilizer and rising food demands from a larger and wealthier human population. These factors point to a global food deficit unless the efficiency and resilience of crop production is increased. The intensification of agriculture has focused on improving production under optimized conditions, with significant agronomic inputs. Furthermore, the intensive cultivation of a limited number of crops has drastically narrowed the number of plant species humans rely on. A new agricultural paradigm is required, reducing dependence on high inputs and increasing crop diversity, yield stability and environmental resilience. Genomics offers unprecedented opportunities to increase crop yield, quality and stability of production through advanced breeding strategies, enhancing the resilience of major crops to climate variability, and increasing the productivity and range of minor crops to diversify the food supply. Here we review the state of the art of genomic-assisted breeding for the most important staples that feed the world, and how to use and adapt such genomic tools to accelerate development of both major and minor crops with desired traits that enhance adaptation to, or mitigate the effects of climate change.

Application of genomics-assisted breeding for generation of climate resilient crops: progress and prospects.

Climate change affects agricultural productivity worldwide. Increased prices of food commodities are the initial indication of drastic edible yield loss, which is expected to increase further due to global warming. This situation has compelled plant scientists to develop climate change-resilient crops, which can withstand broad-spectrum stresses such as drought, heat, cold, salinity, flood, submergence and pests, thus helping to deliver increased productivity. Genomics appears to be a promising tool for deciphering the stress responsiveness of crop species with adaptation traits or in wild relatives toward identifying underlying genes, alleles or quantitative trait loci. Molecular breeding approaches have proven helpful in enhancing the stress adaptation of crop plants, and recent advances in high-throughput sequencing and phenotyping platforms have transformed molecular breeding to genomics-assisted breeding (GAB). In view of this, the present review elaborates the progress and prospects of GAB for improving climate change resilience in crops, which is likely to play an ever increasing role in the effort to ensure global food security.

Genome-wide DArT and SNP scan for QTL associated with resistance to stripe rust (Puccinia striiformis f. sp. tritici) in elite ICARDA wheat (Triticum aestivum L.) germplasm.

KEY MESSAGE: Identified DArT and SNP markers including a first reported QTL on 3AS, validated large effect APR on 3BS. The different genes can be used to incorporate stripe resistance in cultivated varieties. Stripe rust [yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst)] is a serious disease in wheat (Triticum aestivum). This study employed genome-wide association mapping (GWAM) to identify markers linked to stripe rust resistance genes using Diversity Arrays Technology (DArT(®)) and single-nucleotide polymorphism (SNP) Infinium 9K assays in 200 ICARDA wheat genotypes, phenotyped for seedling and adult plant resistance in two sites over two growing seasons in Syria. Only 25.8 % of the genotypes showed resistance at seedling stage while about 33 and 44 % showed moderate resistance and resistance response, respectively. Mixed-linear model adjusted for false discovery rate at p < 0.05 identified 12 DArT and 29 SNP markers on chromosome arms 3AS, 3AL, 1AL, 2AL, 2BS, 2BL, 3BS, 3BL, 5BL, 6AL, and 7DS significantly linked to Pst resistance genes. Of these, the locus on 3AS has not been previously reported to confer resistance to stripe rust in wheat. The QTL on 3AS, 3AL, 1AL, 2AL, and 2BS were effective at seedling and adult plant growth stages while those on 3BS, 3BL, 5BL, 6AL and 7DS were effective at adult plant stage. The 3BS QTL was validated in Cham-6 × Cham-8 recombinant inbred line population; composite interval analysis identified a stripe resistance QTL flanked by the DArT marker, wPt-798970, contributed by Cham-6 parent which accounted for 31.2 % of the phenotypic variation. The DArT marker "wPt-798970" lies 1.6 cM away from the 3BS QTL detected within GWAM. Epistatic interactions were also investigated; only the QTL on 1AL, 3AS and 6AL exhibited interactions with other loci. These results suggest that GWAM can be an effective approach for identifying and improving resistance to stripe rust in wheat.

Genome-wide association mapping for seedling and adult plant resistance to stripe rust in synthetic hexaploid wheat.

Use of genetic diversity from related wild and domesticated species has made a significant contribution to improving wheat productivity. Synthetic hexaploid wheats (SHWs) exhibit natural genetic variation for resistance and/or tolerance to biotic and abiotic stresses. Stripe rust caused by (Puccinia striiformis f. sp. tritici; Pst), is an important disease of wheat worldwide. To characterise loci conferring resistance to stripe rust in SHWs, we conducted a genome-wide association study (GWAS) with a panel of 181 SHWs using the wheat 9 K SNP iSelect array. The SHWs were evaluated for their response to the prevailing races of Pst at the seedling and adult plant stages, the latter in replicated field trials at two sites in Ethiopia in 2011. About 28% of the SHWs exhibited immunity at the seedling stage while 56% and 83% were resistant to Pst at the adult plant stage at Meraro and Arsi Robe, respectively. A total of 27 SNPs in nine genomic regions (1 BS, 2 AS, 2 BL, 3 BL, 3 DL, 5A, 5 BL, 6DS and 7A) were linked with resistance to Pst at the seedling stage, while 38 SNPs on 18 genomic regions were associated with resistance at the adult plant stage. Six genomic regions were commonly detected at both locations using a mixed linear model corrected for population structure, kinship relatedness and adjusted for false discovery rate (FDR). The loci on chromosome regions 1 AS, 3 DL, 6 DS and 7 AL appeared to be novel QTL; our results confirm that resynthesized wheat involving its progenitor species is a rich source of new stripe (yellow) rust resistance that may be useful in choosing SHWs and incorporating diverse yellow rust (YR) resistance loci into locally adapted wheat cultivars.

Optimization of sequence alignment for simple sequence repeat regions.

BACKGROUND: Microsatellites, or simple sequence repeats (SSRs), are tandemly repeated DNA sequences, including tandem copies of specific sequences no longer than six bases, that are distributed in the genome. SSR has been used as a molecular marker because it is easy to detect and is used in a range of applications, including genetic diversity, genome mapping, and marker assisted selection. It is also very mutable because of slipping in the DNA polymerase during DNA replication. This unique mutation increases the insertion/deletion (INDELs) mutation frequency to a high ratio - more than other types of molecular markers such as single nucleotide polymorphism (SNPs).SNPs are more frequent than INDELs. Therefore, all designed algorithms for sequence alignment fit the vast majority of the genomic sequence without considering microsatellite regions, as unique sequences that require special consideration. The old algorithm is limited in its application because there are many overlaps between different repeat units which result in false evolutionary relationships. FINDINGS: To overcome the limitation of the aligning algorithm when dealing with SSR loci, a new algorithm was developed using PERL script with a Tk graphical interface. This program is based on aligning sequences after determining the repeated units first, and the last SSR nucleotides positions. This results in a shifting process according to the inserted repeated unit type.When studying the phylogenic relations before and after applying the new algorithm, many differences in the trees were obtained by increasing the SSR length and complexity. However, less distance between different linage had been observed after applying the new algorithm. CONCLUSIONS: The new algorithm produces better estimates for aligning SSR loci because it reflects more reliable evolutionary relations between different linages. It reduces overlapping during SSR alignment, which results in a more realistic phylogenic relationship.

Characterization of quantitative trait loci controlling genetic variation for preharvest sprouting in synthetic backcross-derived wheat lines.

Aegilops tauschii, the wild relative of wheat, has stronger seed dormancy, a major component of preharvest sprouting resistance (PHSR), than bread wheat. A diploid Ae. tauschii accession (AUS18836) and a tetraploid (Triticum turgidum L. ssp. durum var. Altar84) wheat were used to construct a synthetic wheat (Syn37). The genetic architecture of PHS was investigated in 271 BC(1)F(7) synthetic backcross lines (SBLs) derived from Syn37/2*Janz (resistant/susceptible). The SBLs were evaluated in three environments over 2 years and PHS was assessed by way of three measures: the germination index (GI), which measures grain dormancy, the whole spike assay (SI), which takes into account all spike morphology, and counted visually sprouted seeds out of 200 (VI). Grain color was measured using both Chroma Meter- and NaOH-based approaches. QTL for PHSR and grain color were mapped and their additive and epistatic effects as well as their interactions with environment were estimated by a mixed linear-model approach. Single-locus analysis following composite interval mapping revealed four QTL for GI, two QTL for SI, and four QTL for VI on chromosomes 3DL and 4AL. The locus QPhs.dpiv-3D.1 on chromosome 3DL was tightly linked to the red grain color (RGC) at a distance of 5 cM. The other locus on chromosome 3D, "QPhs.dpiv-3D.2" was independent of RGC locus. Two-locus analysis detected nine QTL with main effects and 18 additive x additive interactions for GI, SI, and VI. Two of the nine main effects QTL and two epistatic QTL showed significant interactions with environments. Both additive and epistatic effects contributed to phenotypic variance in PHSR and the identified markers are potential candidates for marker-assisted selection of favorable alleles at multiple loci. SBLs derived from Ae. tauschii proved to be a promising tool to dissect, introgress, and pyramid different PHSR genes into adapted wheat genetic backgrounds. The enhanced expression of PHS resistance in SBLs enabled us to develop white PHS-resistant wheat germplasm from the red-grained Ae. tauschii accession.

Haplotype diversity of preharvest sprouting QTLs in wheat.

Preharvest sprouting (PHS) is one of the most important factors affecting wheat production worldwide in environments characterized by rainfall and high humidity at harvest. In such environments, the incorporation of seed dormancy of a limited duration is required to minimize losses associated with PHS. A global collection of 28 PHS-resistant and -susceptible wheat germplasm was characterized with microsatellite markers flanking the genomic regions associated with PHS-resistance quantitative trait loci (QTLs), particularly on chromosomes 3D and 4A. The genetic diversity analysis revealed 380 alleles at 54 microsatellite loci, with an average of 7.0 alleles per locus, among the 28 wheat genotypes. Gower's genetic similarity values among all possible pairs of genotypes varied from 0.44 to 0.97, indicating that there is considerable diversity in the PHS germplasm evaluated. Cluster and principal coordinates analysis of genetic similarity estimates differentiated the genotypes into groups, according to their source of PHS resistance. Three major SSR haplotypes were observed on chromosome 4AL, designated RL4137-type allele, Aus1408-type allele, and synthetic-hexaploid-type allele. The RL4137-type allele was prevalent in Canadian cultivars, mostly in cluster 6, followed by the Aus1408-type and its derivatives in clusters 4 and 5. The Syn36 and Syn37 alleles on chromosome 4AL were rare. On chromosome 3DL, the SSRs haplotypes derived from Syn36 and Syn37 were also rare, and proved unique to the Aegilops tauschii - derived synthetic hexaploids. They are therefore likely carrying resistance genes different from those previously reported. Based on genetic relationships, PHS resistance might be improved by selecting parental genotypes from different clusters.

The cre1 and cre3 nematode resistance genes are located at homeologous loci in the wheat genome.

Differential responses in host-nematode pathotype interactions occur in wheat lines carrying different cereal cyst nematode resistance (Cre) genes. Cre1, located on chromosome 2B, confers resistance to most European nematodes and the sole Australian pathotype, while Cre3, present on chromosome 2D, is highly resistant to the Australian pathotype and susceptible to a number of European pathotypes. Genes encoding nucleotide binding site-leucine rich repeat (NBS-LRR) proteins that cosegregate with the Cre3 locus cross hybridize to homologues whose restriction fragment length polymorphism (RFLP) patterns distinguish near-isogenic Cre1 nematode-resistant wheat lines. Genetic mapping showed that the NBS-LRR gene members that distinguished the Cre1 near-isogenic lines were located on chromosome 2BL at a locus, designated Xcsl107, that cosegregates with the Cre1 locus. A haplotype of NBS-LRR genes from the Xcsl107 locus provides a diagnostic marker for the presence of Cre1 nematode resistance in a wide collection of wheat lines and segregating families. Genetic analysis of NBS-LRR haplotypes that cosegregate with Cre1 and Cre3 resistance, together with flanking cDNA markers and other markers from homoeologous group 2 chromosomes, revealed a conserved gene order that suggests Cre1 and Cre3 are homeoloci.